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HYPERSONIC  AREA  RULE 


M.  D.  Ladyzhenskiy  (Moscow) 


Abstract 


Assuming  [2]  that  the  entire  mass  of  gas  is  concentrated  in  the 
infinitely  thin  layer  adjacent  to  the  shock  wave,  report  [ 1] 
formulates  the  hypersonic  area  rule.  According  to  this  rule,  during 
flow  about  thin  nonaxisymmetrical  bluff  bodies  with  egual  values  of 
blunting  resistance  and  identical  laws  of  the  change  in  the  direction 
of  the  flow  on  the  cross-sectional  area,  the  shock  wave  surfaces,  the 
laws  of  the  change  in  pressure,  and,  consequently,  also  the  drag 
acting  on  the  body  coincide.  Here  the  surfaces  of  the  shock  waves  are 
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axially  symmetrical. 

This  study  establishes  the  limits  of  applicability  of  the 
results  in  [ 1 ] and  refines  the  hypersonic  area  rule  by  introducing 
the  entropy  layer. 


1.  Determination  of  the  Limits  of  Applicability  of  the  Results 
of  [1].  As  an  example  of  using  the  hypersonic  area  rule  [1],  wo  will 
construct  a body  equivalent  to  a thin  round  cone,  i.e. , for  which  the 
blunting  resistance  and  the  law  of  the  change  in  the  flow  direction 
of  the  cross-sectional  area  are  identical  tc  those  of  a cone.  The 
cross  section  of  the  body  is  assumed  to  be  shaped  like  an  ellipse 
whose  long  semiaxis  is  equal  to  the  radius  of  the  shock  wave  and 
whose  area  is  equal  to  that  of  the  cross  section  of  the  round  cone 
(Fig.  1).  In  other  words,  the  eccentricity  of  the  ellipse  has  the 
maximum  possible  value  which  still  satisfies  the  requirement 
(condition  (3)  from  [1])  in  each  cross  section,  in  accordance  with 
which  the  body  should  not  go  beyond  the  limits  of  the  space  bounded 
by  the  surface  of  the  shock  wave. 

As  [ 1 ] pointed  out,  the  area  rule  can  be  combined  with  the  law 
of  similarity  during  flow  about  thin  bluff  bodies  f 2 ],  as  a result  of 
which  the  dimensionless  values  which  characterize  the  flow  are 
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determined  by  two  dimensionless  parameters  at  a fixed  value  of  the 
adiabatic  index  * : the  known  similarity  parameter  during  flow  about 
thin  bluff  bodies  K « HooV  *nd  parameter  Kt  » (n  / 2cxs)'l‘ Lx1,  which 
characterizes  the  blunting  effect  of  equal  order  of  aagnitude  to  the 
square  root  of  the  ratio  of  the  resistance  of  the  body  to  the 
blunting  resistance.  Here  r — Sfh/L—  is  the  small  dimensionless  value 
which  characterizes  the  thickness  of  the  body;  S*  is  a certain 
characteristic  cross-sectional  area  of  the  tody;  L is  the  length  of 
the  body;  and  c*,  S are  the  blunting  resistance  coefficient  and  the 
maximum  midsectional  area  of  blunting,  respectively. 

Assuning  that  the  effect  of  blunting  can  te  replaced  by  an 
explosive  effect  at  the  front  point  of  the  body  with  energy  equal  to 
the  blunting  resistance,  blunting  is  insignificant  [1].  Thus, 
blunting  area  S is  introduced  to  express  K,  instead  of  the  blunting 
diameter  [2].  He  will  assume  that  of  the  unperturbed  flow  is  equal 
to  infinity.  Then,  at  fixed  x the  dimensionless  variables  will  depend 
on  parameter  K,  alone. 

Figure  2 (at  *=  1.4}  shows  the  dependences  of  the  value  of  k 
(the  ratio  of  the  long  seaiaxis  of  the  ellipse  to  the  small  one)  and 
(X  - X 0/X o (the  ratio  of  the  resistance  of  the  body  without 
consideration  of  blunting  resistance  to  the  blunting  resistance)  in 
function  Kt  =(«/2c.r5),;,I.tg,a,  in  the  form  of  curves,  where  a is 


the  half-aperture  angle  of  the  round  cone,  lhe  form  of  the  shock  wave 
was  determined  from  the  solution  of  the  problem  of  flow  about  a thir. 
bluff  cone  according  to  [2].  It  is  logical  to  apply  the  area  rule  at 
X/X„  ^ 1.1#  which,  as  Fig.  2 indicates#  corresponds  to  K,  > 0.1.  At 
small  values  of  the  resistance  of  the  body  is  essentially 

» determined  by  the  value  of  blunting  resistance.  At  large  values  of 

Kt,  the  area  rule  loses  its  validity  as  k approaches  one  (more 
precisely  [1],  at  * — 1 ~ (x  — 1) /(x  + l)),  which  occurs  approximately  at  Kj 
= 1.2. 

Thus,  the  range  of  application  of  the  area  rule  lies  within 
* 0. 1 * Kt  4 1.2.  Here  the  ellipse  in  the  cross  section  of  the  body  can 

be  rather  elongated,  differing  from  a circle  (13  k 1.3).  This 
result  can  be  of  practical  interest.  However,  whereas  results  [1] 
were  obtained  with  approximate  assumptions  about  the  concentration  or 
the  entire  mass  of  gas  in  the  infinitely  thin  layer  beyond  the  shock 
wave,  it  is  necessary  to  make  the  area  rule  more  precise. 

S (9 

2.  Refining  the  Area  Rule.  As  before,  we  will  assume  that 
Moa’2>\,  but,  unlike  in  [1],  we  will  not  impose  the  necessary 
condition  t > l.  le  will  introduce  cylindrical  coordinate  system 
xL,  yL,  0 (thn  x-nxis  passes  through  the  front  point  of  the  body  and 
is  directed  along  the  flow).  He  will  use  'vUn,  u,um  to  designate 
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the  velocity  components  in  the  axial,  radial,  and  circumferential 
directions,  respectively:  p?<0^10-~  is  pressure;  pp«>  — is  density; 
is  the  velocity  of  the  unperturbed  flow  directed  along  the  x-axis; 
and  is  its  density.  He  will  characterize  the  blunting  value  by  the 
dimensionless  blunting  diameter  dL,  where  d is  a small  value,  we  will 
write  the  equation  for  the  surface  of  the  body  in  the  form  y = 
r (f)  (x,  «)  . 

He  will  isolate  the  entropy  layer,  i.e.,  the  region  occupied  by 
the  current  lines  passing  through  that  section  cf  the  surface  of  the 
retreating  shock  wave  where  the  angles  of  inclination  formed  by  the 
surface  of  the  shock  wave  and  the  direction  of  the  unperturbed  flow 
are  not  small  (Fig.  3).  Suppose  that  the  equation  of  the  arbitrarily 
introduced  boundary  of  the  entropy  layer  is  y = 6<P  (x,  9),  where  6 is 
small.  Beginning  with  a certain  x = x0  - d,  the  angles  of  inclination 
of  the  boundary  of  the  entropy  layer  to  the  x-axis  will  be  on  the 
order  of  6. 

Below  we  will  give  estimates  of  the  parameters  of  flow  in  the 
entropy  layer  similar  to  those  made  in  [3].  He  will  also  point  out 
that  the  effect  of  the  entropy  layer  on  the  pressure  distribution  on 
a thin  bluff  cone  was  considered  in  [4]. 

He  will  assume  that  the  relationship  / • 


holds  on  the 
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surface  of  the  entropy  layer,  where  a is  the  positive  number  to  be 
deternined. 

As  follows  subsequently,  the  order  of  pressure  is  maintained 
across  the  entropy  layer;  therefore,  with  the  condition  of 

m 

adiaba ticity,  we  can  write  for  density.  Now  we  will  write 

the  continuity  equation  of  the  entropy  layer.  Equating  the  flow  in 
the  entropy  layer  to  the  flow  in  the  jet  current  of  unperturbed  ilow 
through  an  area  equal  to  that  of  the  maximum  midsectional  blunting 
area,  we  will  have 


<2*  — • p u o, 


(1) 


where  * is  the  area  occupied  by  the  entropy  layer  in  cross  section  x 
= const  (t he  shaded  part  in  Pig-  3).  As  follows  from  the  Bernoulli 


equation,  we  will  have  * . 

It  is  obvious  that  estiaate 


since  u * 1 in  the  entropy  layer. 


II 


= J + a, 


(2) 


is  obtained  for  6 in  the  entropy  layer  boundary  equation,  where  s is 
the  cross-sectional  area  of  the  body.  Re  will  stipulate  that  the 
crder  of  nagnitude  of  the  area  of  the  body  does  not  exceed  the  area 
of  the  entropy  layer  S ^ «•  Then,  obviously. 


6 ». 


(3) 


f 
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Since  the  usual  estinate  of  hypersonic  flow  p - 6*  is  valid  for  the 
pressure  on  the  external  boundary  of  the  entropy  layer,  we  will  find 
the  equation  for  determining  a from  it: 


<r~d  . «• 


2x 

* + i 


(4) 


Finally,  we  will  have  the  following  flow  parameters  in  the  entropy 
layer: 


P ~ ax+1.  P ~ ax+l,  6 ~ rfx+1,  a ~ dK+l 


(5) 


with  the  condition  that  the  orders  of  magnitude  of  r and  d are 
related  by  the  relationship  obtained  with  the  condition  s ~ a: 


X ~ dx+l 


(ti 


This  relationship  essentially  conforms  to  the  condition  x-~Vd, 
which  says  that  the  order  of  magnitude  of  the  resistance  of  the  body 
is  comparable  to  the  blunting  resistance  [2], 

We  will  estimate  the  pressure  gradient  in  the  entropy  layer. 

From  the  equations  of  motion  we  have 

dP  _ 


yL  — pfuiLj.,  4-  » 

d»  \ Ox  ' dy  + y gg-J' 

— It-  = — p(uHHL  + v'h’L  4.  *>  d*  \ 

y W do)’ 


O)- 


r 
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Therefore,  J'-J- 


since  all  the  terns  in  the  right  sides  of 
(7)  are  of  the  saae  order  of  smallness.  Whence,  with  consideration  of 
(5)  and  (6),  the  estieate 


»<«+» 


Ap  ~ <fs  — x 


is  valid  for  the  pressure  gradient  in  both  the  radial  and 
circumferential  directions.  Thus,  the  pressure  in  the  entropy  layer 
can  be  considered  to  be  constant  with  the  relative  error 


a 

A P __  rfX  + l __  TX 

p 


(8) 


which  is  somewhat  greater  than  the  relative  error  in  the  theory  of 
small  perturbations  of  the  hypersonic  flow  (equal,  as  we  know,  to 
t*). 


We  will  state  the  flow  problem.  Suppose  that  at  x < x 0 the 
axisy mmetrical  nose  section  of  a body  is  given  with  its  axis  of 


J 
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symmetry  directed  along  the  x-axis,  the  flow  about  which  has  been 
completely  calculated.  As  follows  from  the  above  discussion,  since 
the  entropy  layer  cannot  maintain  the  pressure  gradient  in  ♦■he 
circumferential  direction,  the  pressure  on  cross  section  x = ccnst 
should  be  constant  in  region  x > x0  on  the  outer  edge  of  the  layer. 
It  suffices  for  the  surface  which  bounds  the  entropy  layer  to  be 
axially  symmetrical  in  order  to  satisfy  this  condition.  (Here  its 
equation  can  be  written  as  y = 6 Y ( x)  - ) 


Then  the  flow  outside  the  entropy  layer,  which  is  axi symmetrical 
according  to  condition  at  x < x0»  also  retains  its  axial  symmetry  at 
x > x.  Therefore,  the  condition  of  the  constancy  of  the  pressure  in 
the  circumferential  direction  on  the  external  boundary  of  the  entropy 
layer  will  be  satisfied. 


He  will  single  out  the  equation  which  relates  S and  a at  x > x0, 
for  which  we  will  use  the  continuity  equation,  as  in  [ 3 ].  We  will  use 
the  subscript  0 to  designate  the  values  in  plane  x0.  Isolating  the 
elementary  current  jet  in  the  entropy  layer,  we  will  write  its  flow 
equation: 


PouoJ/odOorfj/o  = puydQdy. 


(9) 


We  will  have  the  following  for  p and  u from  the  adiabatic  and 
Bernoulli  equations,  in  which  the  terms  on  the  crder  of  r*  have  been 


r 
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o = p J-E-Y  ; 4 'iLf-LT*  = 

\ Po  ) l X — 1 Po  V Pd  ) 


We  will  divide  equation  (9)  by  pu  and  integrate  over  the  entire 
entropy  layer  in  cross  section  x = x0  (we  will  use  the  same 
condition:  that  the  coordinates  of  the  selected  current  line  y and  6 
satisfy  relationships  y = y(y0,  90)  , 9 = 9{ y0,  e0)  ) - Obviously,  on 
the  right  side  of  the  equation  we  will  find  area  a occupied  by  the 
entropy  layer  in  cross  section  x.  With  consideration  of  the  axial 
symmetry  of  the  boundary  of  the  entropy  layer,  we  will  have  a = 
v52Y*(x)  - S.  Finally,  the  unknown  relationship  is  written  as: 


(*)  - F (p)  = s (*), 

1 -f.  2 1 _ 2x  Pt 

* — 1 X — 1 Do 

oo 

t + 2 1 _ 2*  Po  (_P_\ 

x — 1 x — 1 po  \po/ 


yorfMyo-  (11> 


Axisy metrical  flov  outside  the  entropy  layer  at  i > x0  can  be 

calculated  by  one  of  the  known  point  methods,  e.g.,  the 
characteristics  method.  Here  the  boundary  of  the  entropy  layer  will 
be  determined  from  equation  (11),  which  plays  a role  in  the  boundary 
condition  which  replaces  the  nonflow  condition,  during  the  solution. 
Whence  it  follows  that  the  flow  is  completely  determined  by  the 
assignment  of  the  law  of  the  change  in  the  cross-sectional  area  of 
the  body  S (X)  at  x > x0.  Here  one  more  obvious  limitation  must  be 
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imposed  on  the  shape  of  the  body:  the  body  must  not  go  beyond  the 
limits  of  the  "entropy  circle"  (see  Fig.  3),  which  can  be  written 
symbolically  as: 

(12) 

Now  we  can  state  the  exact  hypersonic  area  rule  as  follows.  During 
flow  about  thin  bluff  bodies  with  axisymmetrical  nose  sections  which 
coincide  at  a certain  distance  from  the  front  pcint  of  the  body  and 
identical  laws  of  the  change  in  the  cross  section  of  the  remaining 
parts: 

a)  the  flows  outside  the  entropy  layers  are  axisymmetrical;  the 
flow  parameters  at  the  corresponding  points,  the  surfaces  of  the 
shock  waves,  and  the  arbitrarily  introduced  boundaries  of  the  entropy 
layers  coincide; 

b)  the  pressure  in  the  entropy  layers  only  depends  on  x,  and  the 
law  of  the  change  in  pressure  is  the  same  for  the  bodies  in  question; 
as  a result,  the  drags  acting  on  the  bodies  are  equal,  since  drag  X 
is  expressed  as 

i 

* = Xo+L*-PBOUZ,^S'(x)p(x)dx.  (13) 

*. 

where  X0  is  the  drag  of  the  nose  section  of  the  body.  Here  it  is 


■~r-  • 1 
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assumed  that  conditions  (6)  and  (12)  ate  satisfied. 

The  results  obtained  can  easily  be  generalized  to  the  case  of 
flow  with  dissociation.  The  consideration  of  these  phenomena  only 
changes  the  appearance  of  function  F(p)  in  equation  (11). 


3.  Comparison  of  Results.  We  will  compare  the  result  obtained 
with  the  correct  area  proven  in  [1].  The  requirement  of  the 
correspondence  of  the  laws  of  the  change  in  the  cross-sectional  area 
in  the  direction  of  the  x-axis,  as  well  as  the  condition  which  states 
that  the  drag  of  the  body  must  not  exceed  the  order  of  value  of  the 
blunting  resistance,  are  common  to  both  theorems.  The  difference  in 
the  formulation  of  the  theorem  proven  in  §2  is  as  follows: 

a;  the  necessary  condition  m i*  not  imposed  , as  was 
done  in  [ 1 ]; 

b)  instead  of  requiring  the  equivalence  of  the  values  of  the 
blunting  resistance  [1],  a greater  limitation  is  imposed:  the  nose 
sections  of  the  equivalent  bodies,  being  axisy (metrical,  must  have 
the  same  shape; 


c)  instead  of  condition  (3)  in  [1],  which  states  that  the  body 
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cannot  go  beyond  the  limits  of  the  space  bounded  by  the  surface  of 
the  shock  wave,  a greater  limitation  (12)  is  imposed:  the  body  must 
not  go  beyond  the  upper  boundary  of  the  entropy  layer. 

For  this  reason,  we  can  expect  the  values  of  k found  in  §1, 
characterizing  the  difference  in  the  cross  section  of  the  body  from 
that  of  the  equivalent  solid  of  revolution,  to  be  somewhat  higher. 
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